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A THREE DIMENSIONAL MODEL OF THERMOSPHERIC DYNAMICS. 

HI. PLANETARY WAVES 

H. Volland 

Astronomical Institutes 
University of Bonn 
Bonn, Germany 

and 

H. G. Mayr 

Goddard Space Flight Center 
Greenbelt, Md. , USA 

ABSTRACT 

The generation and the propagation of atmospheric waves of planetary 
scale with zonal wave domain numbers n = 1 , 2 , 3 and 4 is calculated within the 
dissipative thermosphere. These waves are generated by solar XUV-heat 
input and by corpuscular heating of the thermosphere. They contain time 
independent zonal components, annual and semiannual components. The 
latitudinal structure of their eigenfunctions and the height variation of this 
structure is determined. It is shown that the relative pressure amplitudes of 
the wave modes approach asymptotic values at high altitudes which are 
proportional to l/n(n+l). Therefore, the thermosphere behaves like a filter 
suppressing the waves with great wave domain numbers n. 

The corresponding wind systems of the various wave modes are deter- 
mined and compared with observed wind data at mesospheric heights as well as 
with density data derived from satellite drag observations. It is shown that 
the predominant annual component with wave number n - 1 is generated by the 
solar radiative heat input which peaks in the summer hemisphere and that the 
semiannual component containing in the mesosphere mainly the wave mode of 
number n = 4 is generated by corpuscular heat input within the auroral ovals. 
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1. INTRODUCTION 

In the first part of this paper, we developed a three dimensional spherical 
model of the thermospherec dynamics (Volland and Mayr, 1971a; referred to as 
paper I). In the second part of this paper (referred to as paper II) we calculated 
the generation and the propagation of three important tidal waves within this 
model thermosphere and compared the results with available experimental 
data. In this third part we want to determine the wave modes of planetary scale 
generated by the corresponding energy modes of solar XUV and corpuscular 
heat inputs. These planetary waves are distinguished by their independence on 
local time (wave domain number m = 0). They can be divided into symmetrical 
terms (zonal wave domain number n is even) which describe the semiannual 
variations and the time independent zonal variations and into antisymmetric 
terms (wave domain number n is odd) which describe the annual variations. 
Terannual and higher order terms appear to be very weak as it follows from 
experimental data and theoretical considerations, therefore, we shall not 
treat them in this paper. 

We shall calculate the equivalent depths and the eigenfunctions of the four 
most important planetary waves with zonal wave domain numbers n = 1,2,3 
and 4. We shall determine the altitude variations in the latitudinal structure of 
these modes and their generation and propagation at thermospheric heights. 
Finally, we shall compare our theoretical result with available experimental 
data. 

The even more important planetary wave mode of wave domain 
number 0, corresponding to an eigenfunction P Q « 1^ will be bypassed 
because this mode does not generate horizontal winds and pressure 
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gradients. This mode has in fact been treated in a one dimensional 
quasistatic model for f « 0 (Jacchia, 1965). The same wave mode with basic 
frequency (*) - H has also been considered in the Harris -Priester-model 
(Harris and Priester, 1962) and for various frequencies by Volland (1969) 
and by Thomas and Ching (1969). 


2. EQUIVALENT DEPTHS AND EIGENFUNCTIONS 

The determination of the equivalent depth and of the eigenfunction of the 
planetary wave modes follows essentially the same line as given in papers I 
and II. The eigenfunction describing the latitudinal structure of the relative 
pressure amplitude is (see (Il/l) ) 


with 



p n f(z) 

P 0 < z > 


0 n *(t?,z) exp 



( 1 ) 


© *(d,z) = ^ f P .(cosd), 

n ' ' n,n' n 1 ' ’ 

n f 

Using the same notation as in (II/l). Here for convenience we 
omitted in (1) the superscript m which is zero for the planetary waves. P r 
are the zonal Legendre polynominals. The planetary waves which we shall 
consider and which are generated by solar XUV -heating and by the average 
corpuscular heat input into the thermosphere have the angular frequency 
a) = 0 (see (1/9) and (1/12)). However, in the case of auroral heating during 
an individual geomagnetic storm , waves qf planetary type can be generated with 
frequencies of the order of owO, where O is the frequency of one solar day 


t 


(Volland and Mayr, 1971b). Therefore, in order to generalize the results 
given below, we shall consider a frequency of 


2fO - Oi + sQ . 

a 


Apparently, it is 


f ~ f - because of 0 a /0 - 1/365, 


where f is the Coriolis parameter introduced in (1/17). 

Contrary to the general behavior of the horizontal winds of the tidal 
waves , the horizontal velocity of the planetary waves must be zero at the 
poles because the horizontal pressure gradient is zero here. Therefore, 
we start from the expressions for the latitudinal and the longitudinal winds 
of (see (II/4)) 


■T = |p exp(2jftlt) 
o F o 


, f P JL 


exp(2jfOt) 


with 


* f =i!^v 5 / 

n 9Ev n* ^n 


rr ^n.n’-l n’-l 


sin# 7 

w V 


0 1 ,P , 

n,n f n' 


Since the dependence on longitude of the planetary waves is zero 

(d/d A =0), the relationship between the coefficients of u and V follows 

from the third equation in (1/16) 

. f _ n* ,f (n'+l) .f 

^n.n* (2n f -l) Vn'-l (2n'+3) ^n.a'+l’ 


( 4 ) 
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Eq. (3) together with (4) indicate that for planetary wave modes with small 
frequency w significant wind velocities can only be expected if the dissipation 
factors in f have significant magnitudes. That is certainly the case at F2 -layer 
heights where due to ion -neutral collisions it is |f y | s 3. However, from wind 
observations at mesospheric heights we have to conclude that even at those 
heights a dissipation mechanism must be effective (see Fig. 1). 

We approximated the wind observation of fig. 1 at Adelaide (latitude: 

35° S) and in Jodrell Bank (latitude : 53° N) by the following series of spherical 
functions 

u= sine j (u x ) o P 1 +(u 3 ) 0 P 3 + ( u 0 P o +u 2 p 2 ) cosf V + ( u 1 P 1 4 V 3 ) coa 2 ° a l J 

V = sinS j (v 0 ) o P 0 + (v 2 ) o P 2 +(v 4 ) q P 4 + (v lPl +v 3 P 3 ) cosO a t (5) 

+ ( v 0 P 0 +v 2 P 2 4v 4 P 4) coa 20 a 1 
Where the coefficients of u and v are related with each other by Eq. (4). 

Using the numbers (wind components in m/sec) 


( U l)o 

= 7.5 ; 

( v o)o 

= 6.2 

( U 3 )o 

= -10. ; 

No 

= 1.8 

u 0 

= -17. ; 

( V 4)o 

= -14.2 

U 1 

• 

o 

II 

v 0 

= 8.3 

u 2 

= 22. ; 

V 1 

= -20.5 

U 3 

= 25. ; 

V 2 

= 43.5 


• 

* 

V 3 

= 33. 


f 


v 


(6) 


= 0.4 j 


= 35.7 
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we plotted our theoretical winds from (5) versus universal time t in fig. 1 
as dashed lines and find satisfying agreement with the observed values 
appart from a prevailing wind component of the order of 10 m/sec in the wind 
data of Adelaine. Note that only the six coefficients u^ are available to fit the 
data and that and ^u^ are responsible only for a shift of the vertical 
scale which is not very well accomplished. Considering the numbers in (6) the 
positive value of the imaginary part of f is remarkable. It can be explained 
readily by eddy viscosity within an atmosphere with a convex vertical profile 
or the longitudinal wind field. (Haurwitz. 1961). Then, 

- fz(" f)~ 2Z vis 0 'V wlth z vi 8 ~ * 0 - 4 < 7 > 

in that height range (see (1/3)). Eddy viscosity of that order and for those 
periods appears to be rather plausible (Kellogg, 1964). Moreover, the mean 
temperature profile is in fact convex between about 85 and 120 km thus sug- 
gesting a similar characteristic for the wind profile there. 

This negative dissipation factor can be effective only below the turbopause 
at about 110km. Above that height, ion-neutral collisions begin to dominate, thus 
leading to positive dissipation factors (see fig. 2)in paper II). Therefore, 
we may expect an effective dissipation mechanism within the lower ther- 
mosphere with a negative value below 110 km, zero at the turbopause and 
positive above that height. 

Having tested the importance of the dissipation factors for the planetary 
waves , we proceed in the determination of their eigenfunctions and equivalent 
depths. We do this in the same manner as outlined in Paper II and look for 
solutions which fulfil exactly the second and third equation in (1/16) and 


s 
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Eq. (II/7). Substituting Eqs. (1), (3) and (4) in the second equation of (1/16), 

f f 

we obtain the following relationship between the coefficients ^ . and 6 • 

n,n n,n 


n*(n' -1 . f 

(2n’ -3)(2n* -1) %,n f -2 


r.2n’ 2 ---2.n'_-l_ . f f 1 J 

L(2n' -l)(2n' +3) VvJ Vn* 


, ( n' +l)(n' +2) f . (2n + 1) t 6 f 

(2n* +3)(2n* +5) ^n,n’ +2 ' ' v=0 n,n’ +1+21* 


(8) 


0 


From (II/7) we obtain 


,,f J . 

. / l y n. n* -1 Mi.n* +l f , x \ *f 

n < n j 2n’ -1 (2n' +3)| TT A n,n’ 6 n,n' _ 


(9) 


with 


!= 1 

n,n’/ 0 


for 


n = n f 
n # n*‘ 


An exact treatment of the eigenfunctions of planetary waves would neces- 
sitate the solution of the infinite set of the equations (8) and (9) provided 
X n n , * 1 for all n* in (9). We shall do this at another place justifying 
there the approximate approach adopted in this paper. 

The coefficients of pressure and horizontal winds in Eqs. (1) and (3) 
have been determined from Eqs. (8), (9) and (4) for the planetary waves with 
wave domain numbers n = 1, 2 , 3 and 4. These coefficients can be written as 

®n,n’ = (°l)n / A n ; * n,n'-l = (°2)n /A n : *n,n’ = / A n 

where the numbers (cj)“ , (cjjj’, (c 3 )“' and * n are given in table 1 . The 
equivalent depths become 


V- 


4f O 


2 2 


n(n+l) f 


( 10 ) 


According to our definition for the normalization of the eigenfunctions 


f f 

0 in (II/6), it is 6 ~ 1. Moreover, for large dissipation factors 

n n,n 


| f u| ’ | f v|» 1 U 18 


6 n.n' « 1/f u f v for n * 

Thus, the eigenfunctions 0 n indeed approach the spherical functions P n 
within the thermosphere: 


0 


0 


n 


n 


( 11 ) 


and the wind functions approach 
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n 


n 




u 


cost? 

zyttj 


U V 


dP 

n 

dd 

dd 


( 12 ) 


above about 200 km altitude. Eq. (12) shows the predominance of the 
latitudinal winds when compared with the longitudinal winds at thermospheric 
heights above 200 km where j z i c i n | > 1 which is due to the dissipation 

mechanisms like viscosity and ion drag. 

Within the nondissipative atmosphere (f u = f = f), the equivalent depths 
in (10) are negative for f < 0. 45 to 0. 7, that is, for periods exceeding 
about one day. The planetary waves are of the evanescent type, there. 

The latitudinal wind is nearly zero for small Coriolis parameters f. 
v/ithin the dissipative thermosphere, the waves behave like quasievanescent 


waves. 




♦ 
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In the next sections we shall discuss the characteristics of the various 
planetary wave modes in detail. 


jT ; 

i: 


p 

r-' 


. 1 1 


2. PRESSURE - AND WIND FUNCTIONS 
2a. Antisymmetric planetary wave modes 

From the XU V-heat inpu* distribution in (1/9) as well as from the exospheric 
temperature distribution in (1/14) follows the existence of the antisymmetric 
planetary wave mode ^ - 0 & /2 O = 1/730^ which at thermospheric 
heights has already degenerated to the zonal spherical function . This 
mode describes solar XUV heat input which peaks in the summer hemisphere and 
it is responsible for the annual shift of the pressure bulge twoard the summer 
hemisphere. From fig. 1 and the numbers in Eq. (6) we have to conclude 

f 

that at least within the lower thermosphere a wave mode 6^ 1 exists too. 

Since the frequency parameter f^ is very small, we can set in the follow- 
ing f - 0 without significant error. Then, the wind components are exactly 
either in phase or in antiphase with respect to the pressure amplitude. The 
phase errors due to our approximation f ® 0 are of the order of 0.5°. 

Fig. 2 gives the pressure function B^ versus co-latitude 9 for the five 
different heights z = 100, 110, 150, 200 and 300 km. Here, we used the 
data of fig. 2 in paper II as the thermospheric parameters and approximated 


L =* = “ JZ, 


k, > 


u v • 'kin ^'kin ®) - 
However* we assumed a dissipation factor = - 0. 4 at 100 in order 


to be consistent with the findings in fig. 1, and we took f » f = 0 at 110 km 
in order to simulate the situation near the turbopause. Fig. 2 indicates 


the usual behavior of the eigenfunctions at thermospheric heights which we 
already observed in the case of the tidal waves in paper II. The eigenfunction 
changes from the Hough -function 0^° at 110 km to the spherical function Pj 
above 200 km altitude. 

In fig. 38 we plotted the wind functions 4^° and 4^° versus co -latitude 
for the height z = 150 km. According to table I, both wind functions do 
not change their latitudinal structure with altitude. Therefore we plotted in 
fig. 3 b the maximum amplitudes of 4^° and 4^° versus altitude. The 
longitudinal wind function 4^° is always negative that is, in antiphase to the 
pressure amplitude. However, the latitudinal wind function 4^° is negative 
below the turbupause and positive (that is , in phase with the pressure) abovo 
the turbopause. At 110 km, the assumed height of the turbopause, the 
latitudinal wind is zero. (If we had taken into account the finite value 
f - f ^ “ 1/730; the latitudinal wind function in fig. 3b would have reached 

_3 

a finite minimum value of 3. 4 x 10 there. ) On the other hand, the 
longitudinal wind function is not affected by the turbopause and reaches a 
maximum near 130 km altitude. Above 160 km, the latitudinal wind exceeds 
the longitudinal wind in magnitude. 

Fig. 4a to 4c give the corresponding pressure and wind functions of the 

- wave i^ode and show essentially the same features as discussed in 
figs. 2 and 3. The wind functions reach maximum values near 150 km. The 
latitudinal wind function 4^° is zero at 110 km. It is in antiphase with the 
longitudinal wind above that height and in phase below that height. 

It should be mentioned that our treatment of the turbulent region below 
the turl pause is subject to some ambiguity. We have set in the numerical 


t 


u 


\ 


f 


12 

calculations of this section f ~ f - 0.4j at 100 km altitude. However, it is by no 
means certain that the curvature of the height profile of the latitudinal wind is 
the same as that of the longitudinal wind. If the sign of f would be opposite to 
the sign c? f , we had to expect according to table I a change of the numerical 
value of A n toward smaller or even toward negative values. Such change can 
influence the amplitudes and phases of pressure and winds. This question 
can be resolved only by a full wave calculation eddy viscosity is inforously 
considered. 

2b. Symmetric planetary wave modes 

From the XUV-heat distribution derived in (1/9), from the distribution of 
corpuscular heating (1/12) and from the distribution of the exospheric tem- 
perature (1/14) we expect the existence of the symmetric planetary wave modes 
0 f 

8 and 8 2 where f 0 = O /Cl = 1/365. Both waves degenerate at thermo- 

spheric heights in',o the spherical function P The wave mode describes the 
time independent zonal heating due to the XUV-energy input which peaks at low 
latitudes and due to corpuscular heating at high latitudes during equinox. The 

f 

wave mode * 8 responsible for the semiannual variation of these heat inputs 
and their corresponding density amplitudes. From (1/12) as well as from the 

f 

findings in fig. 1 we expect also the wave modes 9^ (f = 0 and fg) at thermo- 
spheric heights. 

In figs. 5 and 6 we plotted the usual way the pressure functions 
and the wind functions and <t> n ^ of the two symmetric wave modes 
adopting again the approximation of ~ 0. We observe the usual behavior 
in amplitudes and phases of pressure and winds as discussed in section 2a. 
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As it is apparent from (1/9), (1/12) and (1/13), the most important sym- 
metric planetary wave mode according to our definition is the mode with wave 
domain number n = 0. This mode has no latitudinal structure: 

e o° = p o = 1 • 

and its horizontal winds disappear: 




= 0 . 


It can be treated by an one dimensions ' vertical model as has been done by 
Harris and Priester (1962), Jacchia (1965) and others. 


3. WAVE GENERATION AND PROPAGATION WITHIN THE THERMOSPHERE 
In order to determine the generation and propagation characteristics of the 
planetary wave modes at thermospheric heights , we calculated the normalized 
eigenvalues from Eq. (1/37) and plotted in fig. 7 the attenuation factors 

-1 versus height. They increase with increasing wave domain number n. 

In the case of zero frequency (f-0), all planetary waves behave like pure 
evanescant waves with zero propagation factor or In the case of the annual 
and semiannual wave modes (f - f^ or f^), they become quasi -evanescent with 
propagation factors of the order a n ° ~ 0. 01 which is however of no relevance 
for their propagation characteristics. 

In fig. 8 we calculated from (1/41) and (1/42) the relative amplitudes of 
pressure and vertical wind of the four wave modes generated at thermospheric 
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heights. Here, we adopted the thermospheric parameters at 300 km altitude 

from fig. 2 in paper II and assumed a heat input of 

[-0. 01 z £ 

1 o 


(13) 


0 z < z 

\ o 

As in paper II we want to study the influence of the heat input at ther- 
mospheric heights for the generation of waves. Therefore, the lower boundary 
of the heat input in our model is at z q — 200 km, and the height range z s z q 
shall simulate the situation within the thermosphere above about 200 km. The 
vertical dashed lines in fig. 8 give the asymptotic values of pressure and 
vertical winds according to Eq. (1/45). They have already been reached after 
200 km in the case of wave mode 0^° due to its large attenuation factor 
(see fig. 7). In the case of wave mode 0 1 °, the asymptotic values are still 
twice as great as the amplitudes at z - z q = 200 km. True isobaric layers 

exist for the various wave modes between z-z =8 km (for n = 4) and 

o 

z - z q - 15 km (for n = 1) where pressure and horizontal winds are zero and 
change their phase by 180°. (In the case of finite frequency parameters 
f = f^ or f 2 , a finite minimum value would form there instead, and the phase 
transition would occur within a finite height region. ) Above that isobaric 
layer, the pressure amplitude is greater for smaller wave domain numbers n. 
On the other hand, the amplitudes of vertical velocity are greater for greater n. 
In the case of the vertical winds , the asymptotic value is the same for all 
numbers n as it follows from Eqs. (1/45) and (1/47). The reason for such 
behavior is the following: With increasing wave domain number n, the 
horizontal scale length of the modes decreases. Thus, an increasing amount 


of wave energy is dissipated by viscosity and ion drag which reduces the 
pressure amplitude. On the other hand, the number of the circulation cells 
increase with n (see figs. 9 and 10). In order to drive this wind circulation, 
the vertical velocity must be maintained. Below z q , the height of the lower 
boundary of the heat input, only downgoing waves exist which decay exponent 
ially like 

exp { - (l + 0 n °) (z o - «) / 2H 0 j ; (z<z o ) 


In fig. 9, we plotted the meridional cross-sections of the circulation cells of 
the antisymmetric planetary wave modes and dj* corresponding to the 
heat input of (13) which means heat surplus on the southern hemisphere. 
Longitudinal winds are indicated by the symboles M W" = westerly winds 
(blowing from the west) and "E" = easterly winds (blowing from the east). 

In fig. 10 the same is plotted for the symmetric wave modes 0 2 ° and 0 4 °. 
Here, the heat surplus is at the equator for 0 O ° and at middle latitudes for 



If we relate the wind circulations in figs. 9 and 10 to the annual and 
semiannual winds (f = fj or f^), fig. 9 is valid during the December -solstice 
and corresponds to the heat input 


- J 


2n+l 


cosft t 
a 


Fig. 10 which gives the semiannual variation is valid during the equinox 
and corresponds to the heat input 


| J 2n| 008 20 a l 

Thus, the winds in figs. 9 and 10 change direction after one half period, that is 
after 6 months in fig. 9 and after 3 months in fig. 10. 
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These circulation cells which are identical with thermally driven winds , 
known e. g. from the land -sea -breeze, with rising winds above the maximum 
heat input blowing towards the heat sink above the isobar ic layer and returning 
below that layer should be expected throughout the thermosphere because 
according to figs . 3 and 6 the latitudinal winds of the various wave modes 
although changing their amplitudes with height do not change significantly the 
latitude of their phase transition. The height of the isobaric (or quasi- 
isobaric) layer depends on the real heat input and should lie slightly above 
the maximum heat input which is to be expected generally within the lower 
atmosphere. However, there might be an exception, namely the corpuscular 
heat input which terminates at the base of the thermosphere. If this heat source 
is able to generate a significant wave mode we would expect the isobaric layer 
of the corresponding wind cell near 100 km altitude. This w^uld happen if the 
normalized heat input of the corpuscular heat source exceeds all other 
sources of the same wave domain number n within thermospheric heights. 

It remains to discuss the implications due to the turbopause layer where 
a transitions occurs from eddy viscosity with negative dissipation factors into 
ion drag with positive dissipation factors. In order to maintain the circulation 
cell, we expect a wind field shown schematically in fig. 11. The latitudinal 
wind goes through zero at the turbopause layer without changing its phase. 

On the other hand, pressure and longitudinal wind change their phase abruptly 
although their amplitudes are not affected. In a real atmosphere such transition 
will of course occur within a finite height range. This behavior is quite dif- 
ferent from the situation at the isobaric layer where pressure and both wind 
components become zero and change their phases as indicated in fig. 11. 


4. 


COMPARISON BETWEEN THEORY AND EXPERIMENTAL DATA 


In this section we want to compare the result of the forgoing sections 
with available experimental data. First, we consider the mesospheric wind 
data of fig. 1 and the corresponding coefficients in the series of spherical 
harmonics in Eq. (6). From table I and from Eq. (3), we related the observed 
wind components in (3) to the equivalent relative pressure amplitudes of the 
various wave modes. These numbers are given in table II. Considering these 
numbers , we have to remember our remarks in section 2a about the ambiguity 
concerning the vertical profile of the latitudinal wind. Therefore, the numbers 
in table II should be taken merely as upper limits. Especially, the annual 

* ft 

component 1 may be overestimated due to the sensitive influence of the 
uncertain number = 0. 2 - f u f y £ 0. 36 on the pressure amplitude. 

Nevertheless, the numbers in table II indicate some interesting features. 

In the prevailing wind field (f = 0) the 0 2 ° - component predominates. The 
heat surplus to generate the wind cell of this mode occurs at low latitudes if 
we assume that the height of the mesospheric winds in fig. 1 (92 km) is below 
the isobaric layer of that wave mode. Then, the wind cell of this mode is that 
of fig. 10a. Above the isobaric layer within thermospheric heights the phase 
of the pressure component of this mode is consistent with the corresponding 
T 2 ° - term in the exospheric temperature distribution of (1/14). However 
the amplitude of T 2 ° in (1/14) seems to be underestimated when compared with 
the value suggested from table n. We can at least qualitatively explain this 
discrepancy by assuming that the auroral heat input terminates at about 90 to 
100 km while the solar radiation reaches down to the earth f s surface. There- 
fore the wind system due to auroral heating has its isobaric layer above 92 km 


while the wind system due to XUV heat input has its isobaric layer far below 
that height. Thus in the height range around 90 km both contributions to the 
horizontal wind and pressure fields add together, whereas at thermospheric 
heights they compensate each other due to the reversed signs in the corres- 
ponding heating terms in (1/9) and (1/12). 

The wind direction of the longitudinal wind component at thermospheric 
heights of the Gg^ "mode is consistent with the westerly prevailing winds 
discovered by King-Hele (e.g. , King -He le and Allan, 1966). However the 
amplitude of this wind component estimated from the number in (1/14) and 
from figs. 5c and 8 is 



5 m/sec 


at 200 km altitude and is one order of magnitude smaller than King-Hele's 
wind. Moreover, the longitudinal wind of the - mode disappears at 
the equator which seems to be inconsistent with the findings of King-Hele. 
Therefore, King-Hele's prevailing winds can not be explained by the 
zonal 0g° -mode. 

Even if we assume that Jacchia's (1965) temperature component 
Tg° = - 7° K in (1/14) is an underestimate and that the full contribution of 
the XUV - heat input in (1/9) is effective (neglecting the corresponding 
component of the corpuscular heat input in (1/12)) we arrive at an upper 
limit of the temperature amplitude of the zonal (2,0) mode of 

T 2°U 13 ° K 

as compared with |ti' 1 ^| ~ 120° K of the fundamental diurnal mode. 


That relatively small component of the zonal (2, 0) mode is the result of 
the filtering effect of the thermospheric dynamic system which suppresses 
higher order components of pressure and temperature like l/n(n+l) (see (1/47) 
and (1/49)). This is the explanation for the small observed zonal temperature 
amplitude in (1/14). 

The 0 4 ° -wave mode can be explained as generated by the corresponding 
corpuscular heat input term in (1/12) if the isobaric layer of the wind cell is 
above 92 km, the height of the observed winds in fig. 1. Then, the wind cell 
of this mode hs.3 the form of fig. 10b however with reversed wind directions 
with westerly winds at low latitudes above the isobaric layer. These 
longitudinal winds may be added to the winds of the Q 2 ° -mode. However, 
at high latitudes the winds of this mode are easterly. The maximum westerly 
wind of the $ 4 ° -mode at 200 km altitude is estimated from (1/12) and figs. 6c 
and 8 as 

(v) 4 °| * 1 m/sec 


which again is too small to explain the King-Hele-wind. The relative pressure 
amplitude of this mode at thermospheric heights is expected from the cor- 
puscular heat input in (1/12) and already transformed into Jacchia's (1965) 
exospheric temperature amplitude is 


s 2° K 


(T q * 1000° K) 


which appears to be not very significant. 

£ 

We now turn to the annual components in table n. Here, the 6^ 1 - com- 
ponent predominates although the number in table II may be overestimated as 
already mentioned. This component is consistent with the corresponding 


XUV-heat input term in (1/9) and with the exospheric temperature term in 

(1/14) if the i8obaric layer of this mode is well below 92 km, the height of 

the wind in fig. 1. Then, the heat surplus occurs on the summer hemisphere 

and can readily be explained as resulting from the solar radiative heat input 

which peaks in the summer hemisphere during solstice. 

From (1/9) and from the result of section 3 we estimate a pressure 

component related to the pressure of the fundamental symmetric diurnal 
1 1/2 

mode 8 j ’ (see paper II, section 3) of 


Pl f ‘ 


3 l' 

1,1/2 


.1,1/2 

Pi 


J 1 


which agrees completely in amplitude and phase with the corresponding number 
of the exospheric temperature distribution of (1/14): *) 


ft 


xJn 


0. 048 

0. 12 


0.4 


*) Note that Jacchia's exospheric temperature (Jacchia 1965) represents 

merely the density amplitude above about 200 km. A ratio of T n m /T Q = 0. 1 

corresponds to a relative density amplitude of p m /p ~ 0. 4 for T = 1000° K 

no o 

and above 200 km. Moreover, the density amplitude is proportional to the 

pressure amplitude, and it is p^ m /p^ ~ 0. 85p„ m /b within thermospheric 

no n o 

heights above about 200 km. 
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The wind cell of fig. 9a is just the wind circulation proposed by Kellog 
(1961) and by Johnson (1964) to explain the F2 -layer winter anomaly and the 
winter helium bulge. Moreover, from our result in section 2, from fig. 11 
and from (1/42) it follows that within the turbulent region with convex wind 
and temperature profiles the pressure and temperature are out of phase with 
the latitudinal wind. Therefore, in the height region near 90 km we expect a 
cold summer pole and a warm winter pole which is consistent with experimental 
data (Cole and Kantor, 1963). The amplitude of this observed temperature 
wave is of the order of Tj 1 /T ~ 0. 25 at 90 km altitude which is quite con- 
sistent with the number given in table II. 

We estimate the maxima latitudinal wind at 200 km from the number in 
(1/9) and from figs. 3b and 8 as 

f I 

(u)j 1 ~ 25 m/sec 

which is a reasonable value to explain the F2 -winter anomaly. 

f f 

The 1 -component in table n has the opposite sign of the 0^ 1 -component. 

Therefore provided the heat input maximum occurs below 90 km, its corres- 
ponding wind cell has the form of fig. 9b, however, with opposite wind direc- 
tion. That is equivalent to a heat surplus on the winter pole. The corres - 
ponding XUV-heat input term is very small {in our approximation in (1/9) it is 
exactly zero). The origin of this component seems to be obscure, therefore. 

It may perhaps be related to the real heat distribution within the lower 
atmosphere where such term would damp the pure cos0 - dependence of the 
annual heat variation as the result of some horizontal heat transfer from the 
summer to the winter hemisphere. 
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Next, we consider the semiannual components in table 11. The latitudinaly 
depending part of the semiannual effect at thermospheric heights is according 

f 

to table II predominantly included in the 0^ 3 -mode. Its phase in table II is 
consistent with the corresponding corpuscular heating term in (1/12) if we 
assume that the isobar ic layer of that mode is above 92 km, the heights of the 
winds fig. 1 . This behavior , therefore , can be considered as typical for 
wind systems driven by the corpuscular heat input. The wind system of this 
mode is that of fig. 10b, however, with reversed wind direction. 

We estimate from the number in (1/12) and from the results of the 
foregoing sections a pressure amplitude of this mode above 200 km already 
transformed into the exospheric temperature amplitude of 



which is completely insignificant at thermospheric heights. 

f 

The small © 2 2 component at 92 km (see table II) does not fit into the 
general picture outlined above because it should be relatively large as 

f 

compared with the 3 component due to the additive effects of solar radiation 
and auroral heating at this height range (see the remarks in connection with 
the zonal © 2 ° component). However at thermospheric heights we expect rather 
small amplitudes of the O^ 2 component with an upper limit of the temperature 
amplitude according to (1/9) of 

3° K (T o = 1000° K), 

which is not inconsistant with Jacchia's analysis in (1/9). 
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Finally, we want to estimate the amplitudes of the zonal and the semi- 
annual (0,0) modes generated by solar XUV and auroral heating. We do this 
by adopting the half empirical formula (1/49) which is valid for n - 0 and 
which gives rather good approximated values at thermospheric heights as has 
been tested by full wave calculations. Since the fundamental diurnal mode 


1 1/2 

6^’ has the well established temperature amplitude of 


, 1 , 1/2 


~ 120° K 


during moderate solar activity (T q = 1000° K) we shall use this value as a 
normalizing factor. A second normalizing factor will be adopted from the 
upper limit of the corpuscular heat input in (1/15). Then we determine from 
(1/9), (1/12), (1/15) and (1/49) a zero component of the exospheric temperature 
of 



2.4 x 1. 17 x 120 ,_o „ 

' T4TT83 460 K 


(14) 


as the mean temperature at thermospheric heights due to the contributions of so- 
lar XUV and auroral heating. This temperature amplitude is smaller by a factor 
of about two when compared with the observed mean exospheric temperature 
of T q - 1000° K in (1/14). The number in (14) is qu!*e consistent with the 
finding of Jacchia (1965) that the observed mean exospheric temperature 
extra - polated to zero solar activity is about 500° K. 

We therefore need an additional heat source to generate the residual 500° K 
temperature. Several heat sources have been proposed. Hines (1965) assumed 
heating due to wave dissipation of the whole spectrum of gravity waves within 
thermospheric heights. Volland (1969) assumed wave dissipation of the funda- 
mental diurnal mode. Lindzen and Blake (1970) proposed wave dissipation of 
the semidiurnal modes. Since the semidiurnal waves are propagation modes 
within the lower atmosphere with relatively large amplitudes near 100 km 
height, Lindzen and Blake's (1970) assumption is probably the most likely one. 


ii. 


; 



)( 
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The mean heat input due to that dissipated energy in (1/13) has the 
amplitude 

Q dis" L4Q XUV < 15 > 

during moderate solar activity (T Q = 1000° K). 

Turning to the semiannually varrying (0, 0) component in (1/12) we arrive 
at an amplitude of the corpuscular heat input of 

T 0 f J | s 13° K. (16) 

From Jacchia's analysis (1/14) we find an amplitude of 

T 0 f2 | = 40° . (17) 


Therefore an additional heat source must be effective to generate the tem- 
perature difference which is again most likely due to wave energy dissipation of 
tidal waves (Volland, 1969). We estimate its energy amplitude from the 
numbers in (14) to (17) using the already mentioned assumptions as 



diss 

which is just the value proposed in (1/13). 


0.05 



5. CONCLUSION 

We considered in this paper planetary waves at thermospherics heights 
generated by solar XUV-heat input and by corpuscular heating. We determined 
the eigenfunctions of the four most important planetary waves with zonal wave 
domain numbers n = 1,2,3 and 4 and showed the height dependence of the 
latitudinal structure of these eigenfunctions. Due to dissipation forces like 
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viscosity anol ion drag, the eigenfunctions of the wave modes become 
essencially identical with the zonal spherical harmonics P n (cos0 ) 
above about 200 km. 

We calculated the generation and the propagation of the various wave 
modes within the dissipative thermosphere. It has been shown that the 
relative pressure amplitudes and the horizontal winds of the wave modes 
approach asymptotic values at high altitudes which are proportional to 
(l/n(n+l). Therefore, the thermosphere behaves like a filter for the pres- 
sure amplitudes and the horizontal winds suppressing waves with large zonal 
wave domain numbers n. However, the amplitudes of the vertical winds have 
the same asymptotic value for all wave modes. 

The corresponding wind systems of the various wave modes have been 
determined and compared with observed wind data at mesospheric heights 
as well as with density data derived from satellite drag observations. The 
following conclusions have been drawn: 

1. The prevailing winds (zero period) containing the components 
n = 2 and 4 are rather strong at mesospheric heights, however small at 
thermospheric heights. The 0 2 ° -component has its origin in solar radiation 
and in auroral heating. Both heating components which are in antiphase 
contribute additionally to the winds and pressure fields near 90 km because 
the auroral heat input which terminates near 100 km altitude generates a 
wind cell with an isobaric layer above 90 km. On the other hand, the solar 
radiation generates a wind cell with an isobaric layer below that height. At 
thermospheric heights, however, the winds and pressure fields of both heat 


sources compensate each other. The 6^ -component can be explained as 
generated mainly by auroral heating. 

2. The annual winds with components n = 1 and *3 are strong at 

f 

mesospheric heights. The 6^ 1 -component is also significant at thermo- 
spheric heights and can be explained by the solar radiative heat input which 
peaks in the summer hemisphere during solstice. Due to eddy heat conduc- 
tivity at mesospheric heights below the turbopause , a cold summer pole near 
90 km altitude is expected from our theory wnich is consistent with observa- 

f 

tions (Cole and Kantor, 1963). The 6^ 1 -component is strong at 90 km, how- 
ever supressed already at thermospheric heights due to the filtering effect of 
the thermosphere. Its origin is obscure. 

f 

3. The semiannual wind component 9 2 is insignificant at 90 km as well 

2 

as at thermospheric heights. This is e^lained by the compensating effect in 

f 

the corresponding components of XUV and corpuscular heating. The 6^ 2 -com 
ponent predominates the latitudinal structure of the mesospheric winds near 90 
km. It is explained as originating from the corresponding heating term in the 
auroral heat input. 

4. The mean exospheric temperature can be considered as the tempera- 
ture amplitude of the zonal wave mode with wave domain number n = 0 and 
zero frequency (f - 0). It is generated by solar XUV heat input and to a 
significant amount by the dissipated wave energy due to tidal waves from the 
lower atmosphere. The semiannually varrying zonal wave mode oj 2 has 

its origin mainly within the lower atmosphere and results from wave dis- 
sipation of planetary waves from below. 
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Table II 

t 

Relative pressure amplitudes p R /p Q of the various planetary wave 

f 

modes 6 n deduced from the mesospheric winds observed at 92 km 
altitude. (f=0 : prevailing component ; f ~ f^ : annual component ; 
f = f. : semiannual component) 


x 


1 

2 

3 

4 


0 

-0.082 

0 

0.032 


0. 272 
0 

•0. 098 
0 


0 

0.005 

0 

-0. 081 
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Figure 1. 

Figure la. 
Figure lb. 
Figure 2. 

Figure 3a. 
Figure 3b. 
Figure 4. 
Figure 5. 
Figure 6. 
Figure 7. 


FIGURE LIST 

Mesospheric winds at 92 km height observed at Adelaide (35° S) and 
Jodrell Bank (53° N) versus time of year (solid lines) after 
Kochanski (1963) and theoretical winds (dashed lines). 

Meridional wind (N'ly ~ winds blowing from the north). 

Longitudinal winds (W'ly = winds blowing from the west). 

Pressure function of the antisymmetric planetary wave mode 0^° 
versus co -latitude & for five different heights within the 
thermosphere. 

Latitudinal wind function 4^ and longitudinal wind function 4>^ 
versus co -latitude d at 150 km altitude. 

Maximum wind amplitudes of the 0^° -mode versus altitude. (Solid 
line: positive; dashed lines: negative). 

Pressure and wind functions of the antisymmetric planetary wave 
mode 6^ versus co-latitude. 

Pressure and wind functions of the symmetric planetary wave mode 
0 2 ° versus co -latitude. 

Pressure and wind functions of the symmetric planetary wave 

mode 0.^ versus co -latitude. 

4 

Attenuation factors -1 is the imaginary part of the 
eigenvalue of wave modes with wave domain number n. ) 
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Figure 8. 


Figure 9. 


Figure 10. 


Figure 11. 


Relative wave amplitudes of pressure p and vertical wind w of 
the planetary wave modes versus altitude within the dissipative 
thermosphere (z q ~200 km). The generating heat input is in phase 
with the pressure above the isobar ic layer (solid lines) and is in 
antiphase with the pressure below that height (dashed lines). It 
terminates at z q . Below that height, only downgoing free internal 
waves exist. 

Meridional cross-section of the wind cell of the antisymmetric 

planetary wave modes with annual period fj during December - 

f f 

solstice. (Fig. 9a: 8^ 1 -mode; fig. 9b: $ 3 1 -mode.) 
Longitudinal winds are indicated by W - westerly and E = easterly. 
Meridional cross-section of the wind cell of the symmetric 
planetary wave modes with semiannual period f 2 during equinox 
conditions. (Fig. 10a: 8 -mode; fig. 10b: 8* 2 -mode.) 
Longitudinal winds are indicated by W = westerly and E = easterly. 
Schematic picture indicating amplitude and phase relationship 
between the wind components near the turbopause and near the 
isobaric layer. Equator is on the left, poles are on the right. 
Longitudinal winds (E = winds from the east; W = winds from 
the west) in brackets yield for the southern hemisphere. 






PRESSURE G 







HORIZONTAL WIND 
OF 0°- MODE 

MAXIMUM WIND 




Figure 3a. Latitudinal wind function 0 and longitudinal wind function ‘i»,° versus 
co-latitude # at 150 km altitude. 

FigureSb. Maximum wind amplitudes of the 0 , 0 mode versus altitude. (Solid 
line positive, dashed lines negative). 
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F igure 5. Pressure and wind functions of the symmetric planetary wave mode 0 2 ° versus co-latitude. 
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Figure 6. Pressure and wind functions of the symmetric planetary wave mode 0 4 ° versus co-latitude. 
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Figure 7. Attenuation factors £f n ° -1 (0 n ° is the imaginary part of the 
eigenvalue of wave modes with wave domain number n.) 












CIRCULATION CELLS OF 
ANTISYMMETRIC PLANETARY MODES 
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